Shot Noise of a Tunnel Junction Displacement Detector 
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We study quantum-mechanically the frequency-dependent current noise of a tunnel-junction cou- 
pled to a nanomechanical oscillator. The cases of both DC and AC voltage bias are considered, 
as are the effects of intrinsic oscillator damping. The dynamics of the oscillator can lead to large 
signatures in the shot noise, even if the oscillator-tunnel junction coupling is too weak to yield an 
appreciable signature in the average current. Moreover, the modification of the shot noise by the 
oscillator cannot be fully explained by a simple classical picture of a fluctuating conductance. 
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Spurred primarily by experiments in solid-state qubit 
systems, there has recently been considerable inter- 
est in understanding the noise properties of mesoscopic 
systems used as detectorai^^^. Many new re- 
sults have emerged, including an understanding of the 
connection between noise, back-action dephasing and 
information™^, and of the influence of coherent qubit os- 
cillations on the output noise of a detector—. Not surpris- 
ingly, similar concerns arise in the study of nanomechan- 
ical oscillators. Recent experiments using single-electron 
transistors (SETs) have demonstrated displacement de- 
tection of such oscillators with a precision close to the 
maximum allowed by quantum mechanics^. Given the 
interest in these systems, it is important to gain a better 
understanding of how a mesoscopic detector influences 
the behaviour of an oscillator, and vice- versa. Several 
works have addressed various aspects of this problem. In 
particular, it has been shown that an out-of-equilibrium 
detector can serve as an effective environment for the 
oscillator, providing both a damping coefficient and an 
effective temperatureiiiiSii 3 .. 

In the present work, we turn our attention to the finite- 
frequency output noise of a mesoscopic displacement de- 
tector, where one expects to see signatures of the time- 
dependent fluctuations of the oscillator. A completely 
classical study of the current noise of a DC-biased SET 
displacement detector was presented recently in Refs. Il ll 
andE In contrast, we consider a generic tunnel-junction 
or quantum point-contact (QPC) detector, in which the 
tunnelling strength depends on the position of the os- 
cillator, and calculate quantum mechanically the finite 
frequency current noise. Such a system could be real- 
ized by using an STM setup where one electrode is free 
to vibrat oV . We treat both DC and AC voltage bias; 
the latter is of particular interest, as in experiment, it is 
common to imbed the detector in a resonant tank circuit 
for impedance-matching purposes, and then probe its AC 
response. We find that even for a detector-oscillator cou- 
pling so weak that there is little signature of the oscil- 
lator in the average current, there can nonetheless be 
a strong signature in the finite-frequency current noise. 
We moreover find that the oscillator contribution to the 
noise cannot be simply explained by a classical model of a 



detector conductance which fluctuates with the oscillator 
position- there are additional quantum corrections which 
suppress the contribution of zero point fluctuations. We 
show that these quantum corrections result from corre- 
lations between the detector's random back-action force 
and intrinsic noise. Finally, in the AC-biased case, we 
find that the oscillator experiences a time- dependent tem- 
perature, which has a direct influence on the detector's 
current noise. 

Model- Considering the simplest case where the tunnel- 
matrix element depends linearly on the oscillator dis- 
placement x, the tunnel junction detector is described 
by: 
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Here, CL,k (cfl.fe) destroys an electron state in the left 
(right) electrode, A is the conduction-electron density of 
states, rh denotes the number of tunnelled electrons, and 
the operator Y< augments rh by one. r\ parameterizes 
the sensitivity of the transmission phase to x, and will 
in general be non-zeroi&. We consider both the cases 
of a pure DC voltage, V(t) = V and a pure AC volt- 
age, V(t) = V cosvt. Note that the tunnelling Hamil- 
tonian itself acts as a random back-action force F on 
the oscillator; this corresponds to random momentum 
shifts imparted to the oscillator by tunnelling electrons 9 . 
We will describe our system by a reduced density matrix 
p{m\ x, x'\ t) = (x\p(m; t)\x') which tracks the state of 
the oscillator and m, the number of electrons which have 
tunnelled through the junction. As there is no supercon- 
ductivity, p is diagonal in m. In general, the evolution of 
p will be given by a Dyson-type equation: 



dt 



p(m,t) = -UH ,p(m;t)}+ / dt'J2 
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Here, Uq is the evolution operator corresponding to the 
unperturbed (zero-tunnelling) Hamiltonian, and we have 
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FIG. 1: Diagrams for the a) scattering in, and b) scattering 
out terms in the self energy £(£— t'). The solid lines represent 
the forward and backwards Keldysh contours; the dashed lines 
are conduction electron propagators. The solid black vertices 
correspond to to + t'x. Note that an x operator appearing in 
a vertex at time t' will evolve during the tunnelling event. 

written the self-energy E as a super-operator (i.e. an 
operator acting on the space of density matrices). 

We will consider the simplest case of weak tunnelling, 
and keep only self-energy terms which are lowest order 
in the tunnelling. E is only non-vanishing if m! = m 
or m! = m ± 1; these two types of contributions corre- 
spond to "scattering out" and "scattering in" terms in a 
kinetic equation, and are given by the diagrams shown 
in Fig. 1. These diagrams correspond to standard tun- 
nelling bubbles^, the only difference being that the tun- 



nelling vertices can contain an x operator. If x appears 
at the t' end of a graph for E (£,£'), x will evolve during 
the duration of the tunnelling event. As a result, the 
self energy E has terms involving p, and the final form 
of E we obtain does not correspond to the oscillator-free 
case with x dependent rates. We also include perturba- 
tively the effects of a high-temperature Ohmic heat bath 
(ksTbath 3> h£l, with Q being the oscillator frequency) 
on the oscillator using a Caldeira-Leggett description 18 
and the lowest-order Born diagrams in the self-energy 
(i.e. same diagrams as in Fig. 1, with tunnelling bub- 
bles replaced by environmental boson lines). 

Finally, we specialize to the case where the voltage 
V is much larger than Ml/e. For weak tunnelling, 
eV 3> Ml, and small AC frequency v, it is then possi- 
ble to make a Markov approximation in Eq. : U Q (t — 
t')p{m']t')ul{t — t') — > p(m';t). We are assuming that 
over the short timescales relevant to tunnelling, one can 
describe the dynamics of the density matrix by its zero- 
tunnelling evolution. Fourier-transforming in the m in- 
dex, p(k;t) = J2m=~oo e lkm p(m; t), Eq. becomes: 
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where 70 is the intrinsic damping coefficient associated 
with the equilibrium bath, Dq = 2M^fkBT(, a th is the 
corresponding diffusion constant, and a = +(— ) labels 
contributions from forward (backwards) tunnelling. The 
detector-dependent diffusion constant D(t) = J2a-D<j(t) 
and damping coefficient j(t) = X^o-7o-(i) are given by: 

- M3\ m -™- ,a > 

D ° {t) = (r<r( *' m )+ T "( t >- hQ )) ( 4 ) 

while F(t,rf) = sin r/ (^^j 2o~D a (t)/h is the average 

back-action force exerted on the oscillator. T±(t,E) are 
the t 1 = finite temperature forward and backwards 
inelastic tunnelling rates involving an absorbed energy 
E; these rates are time-independent in the case of a DC 
voltage. Note that we have neglected self-energy terms 
which rcnormalize the oscillator Hamiltonian; these are 
unimportant in the weak-tunnelling limit we consider. 

Eq. @ yields a compact description of the cou- 
pled detector-oscillator system; it is a generalization of 
an equation first derived (via an alternate approach) 



by Mozyrsky et a/. 11 to an arbitrary detector in the 
tunnelling regime, including the possibility of an in- 
dependent tunnelling phased, a nonlinear junction I- 
V, a time-dependent bias voltage, and intrinsic oscilla- 
tor damping. Taking k = yields the equation for the 
reduced-density matrix of the oscillator, and (c.f. Ref. 
ITlh has the Caldeira-Leggett form for a forced, damped 
oscillator in the high-temperature regime^. In what fol- 
lows, we focus for simplicity on the case of T = in the 
tunnel junction, and on 77 = 0, which ensures F = 0i&; a 
non-zero F does not significantly change our results. 



Shot Noise- Eq. (J2J can in principle be used to calculate 
the full counting statistics of tunnelled charge as a func- 
tion of time. By focusing solely on the time-dependence 
of the reduced second moment ((m 2 (t))) (i.e. variance), 
it is possible to calculate the symmetrized frequency- 
dependent current noise using the MacDonald formula^. 
In the case of an AC bias voltage, the noise will be a func- 
tion of two times. We focus on the part of the noise that 
is independent of the average time co-ordinate, a quan- 
tity which is directly accessible in experiment. It is given 
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by a modified version of the MacDonald formula: 



Si(u>) = 2e 2 ui 



dt sin Lot 



2n 



d t ((m 2 (t^))) (5) 



where </> is the initial phase of the AC voltage. 

DC Bias- For a DC biased normal-metal junction at 
T = 0, the tunnelling rates are given by hr a (t,E) — 
{T ) 2 (aeV + E)Q(aeV + E). Eqs. ©-0 yield 7 = 
ftr' 2 /(47rM) and k B T eff = eV/2 11 , We find from Eqs. 
and © that the current noise may be written as 
Si(uj) = 2e(I) + ASi, where the first term corresponds 
to purely Poissonian statistics, and the second term is a 
correction arising from correlations between the motion 
of the oscillator and the number of tunnelled electrons: 



may nonetheless be a large peak in the noise if eV / (h'jtot) 
is large. The upper bound in Eq. JSJ) corresponds to the 
optimal scenario, where there is no intrinsic (detector- 
independent) damping, and a> 1. The maximum S/N 
is determined by eV and the sensitivity t'/to, and can 
be arbitrarily large. Due to the dependence on 7, we find 
the surprising result that the maximum S/N is inversely 
proportional to the detector sensitivity t' /to. Note the 
marked difference from experiments attempting to detect 
coherent qubit oscillations in the detector current noised, 
where back-action effects limit the S/N to a maximum 
of 4. 

We turn now to the second term in Eq. 0, which is 
a lower-order in V quantum correction to the classical 
result. It would appear to cause AS/|i to vanish in the 



AS! (to) 



Ae A V 



dt sin Lot 
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(2r r') ((x(t) ■ m(t))) + (r'f ((x 2 (t) ■ m(t)» 



Physically, the covariances appearing above arise from 
the x-dependence of the tunnelling probability- if m(t) 
is larger than average, then it is likely that x(t) and x 2 (t) 
are also larger than average. These covariances can be 
calculated directly from Eq. © , and obey simple classical 
equations corresponding to a forced, damped harmonic 
oscillator. Consider first the contribution from {{x ■ m)) 
in Eq. ©, which is leading order in t'. In calculating this 
covariance, one finds that the tunnel junction provides an 
effective driving force; we find a contribution: 



AS7(w) 
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= — (2r r) 



S x (u>) (7) 



where S x (w) = 8 7 ft 2 (x 2 )/((w 2 - fl 2 ) 2 + 4 7 2 w 2 )- 1 is the 
spectral density of oscillator x fluctuations obtained from 
Eq. (J2J), and (Axo) 2 = h/(2MQ) is the zero-point un- 
certainty in the oscillator position. The first term in 
Eq. Q is exactly the answer expected (to lowest order 
in t') from a simple picture of a classically fluctuating 
junction conductance (i.e. ASi(lu) = V 2 Sg{uj), where 
Sg(u)) is the spectral density of conductance fluctuations, 
and is in turn determined by S x (lo)). Equivalently, if 
we think of our junction as an x-to-I amplifier having 
a gain A = 2eVr r'/h, this first term corresponds to 
simply amplifying up the fluctuations of the oscillator: 
ASi = X 2 S X . Eq. Q yields a peak in 5/(w) at lo = O; 
keeping only the leading term in V, the ratio of the peak- 
height to the background Poissonian noise (i.e. the S/N 
ratio) is: 



(8) 

where a 2 = t' 2 (x 2 )/tq, jtot = 7o + 7- Note that if a is 
small, there will be no sizeable signature of the oscillator 
in the average current (i.e. S{I)/(I)q ~ a 2 ), but there 



limit eV 



m/2, {x 2 



(Axq) , i.e. it suppresses 



a zero-point contribution to A«S7|i. (Of course, we can- 
not rigorously take this limit, as Eq. is strictly only 
valid for eV ^> Ml). A similar result was found for the 
average current (7) in Ref. 0, where a similar offset 
term could be traced to the inherent asymmetry between 
events in which energy is absorbed from the oscillator, 
versus those in which it is emitted to the oscillator. In 
the present case, the quantum correction to the noise in 
Eq. |J7J| can be given a classical interpretation- it arises 
from correlations between the intrinsic shot noise of the 
detector, and the back-action force F acting on the os- 
cillator. If there are such correlations, we would expect 
classically the current noise to have the form: 

AS7(w) = X 2 S x (oj) + 2ARe [g(-uj)S IF (uj)} (9) 

where Sif(uj) is the symmetrized cross-correlator be- 
tween the junction current and back-action force noise, 
and g(oj) is the oscillator response function. Note that 
the second term above is oc V, while the first is cx V 2 . 
It is well known that inversion symmetry forces Re Sif 
to vanish; this is what allows a QPC detector to reach 
the quantum limit for measuring a qubils^. However, 
at finite lo, Im Sif is non-zero. Consequently, the sec- 
ond term above is non-zero; a direct perturbative cal- 
culation (assuming a thermal state for the oscillator) 
shows that this term corresponds to the second term in 
Eq. |0| . Thus, we see that quantum corrections to the 
noise, which suppress zero-point contributions, can be as- 
sociated with classical out-of-phase correlations between 
the random back-action force and the intrinsic detector 
output noise. 

Finally, we return to Eq. © and examine the contri- 
bution from ((a: 2 • m)), a term which is higher-order in 



One finds: 



e 3 V , 4 eV n (Ax ) 2 



[ASj(w)] 2 = —{r'f -r 



h v ' ' \h 2ir (x 2 ) I 

x J ^S x {J)S x {uo-J) (10) 

Again, the first term above agrees with the expecta- 
tion for a classically fluctuating junction conductance; 
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FIG. 2: Oscillator contribution to Sf, the frequency- 
independent part of the shot noise (i.e. second term in 
Eq. 11310 . versus the AC voltage frequency v, for Sl/7 = 
50, 70 = 0, and eV 3> hv, hSl. The maximum suppression 
of this term at v = £7 (over its v — > 00 limit) is by 8/9. 
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FIG. 3: Full shot noise Si (to) for an AC bias voltage of fre- 
quency v — 100£7, including the effects of correlations between 
xCt) and m(t) (c.f Eq. JT1J). We have chosen eV = 100?ifi, 
a 2 = 1, (70 + 7) = £1/20, and to = 0.1. The y-axis is scaled 
by the value of the ^-independent part of the noise. 



it yields peaks in Si at u> — and u = 20. The second 
term is a quantum correction, completely analogous to 
that found for ASVJi- 

AC Bias- We now consider an AC bias voltage 
V(t) = Vcos(vt), where eV 3> hv,h£l. In the limit 
of small v, it is possible to derive a simple expres- 
sion for the time-dependent tunnelling rates^S. Defin- 
ing hf(E) = {t ) 2 E ■ O(E), we have T a {t,E) = 

En=o i 1 - <W/2) v n T { ?\E) cos nut, with: 



p(n) 



■ 



7T ^ / 

— cos(n(?)r ( eFcos^ 

7T V 



„ nhv\ , 
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Using Eqs. ©-101, we find that the damping coefficient 
7 of the oscillator is time-independent and identical to 
that in the DC case, whereas the diffusion constant is 
time-dependent and contains higher harmonics of the AC 
frequency v. Writing D(t) = 2MjksT e ff(t), we have to 
a good approximation: 



k B T eff (t) 



,Y rg" ( 2(-l)« 
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The small but finite photon frequency f prevents higher 
harmonics from contributing to T e ff] without it, we 
would have simply ksT e ff[t) — V\ cosi4|/2, which tends 
to zero twice each period. With the finite cut-off in- 
cluded, the minima of kBT e ff(t) are ~ hv. The time- 
dependence of T e f f (t) implies that the position variance 
(x 2 (t)) of the oscillator will be time-dependent and in 
phase with the AC voltage; as we show, this has a direct 



influence on the noise and the average current. For the 
latter quantity, we find: 



(Ht)) 



,2y 



cos(^) (T 2 + (r')V(*)>)-AI(t) 



where the quantum correction is approximately A/(t) ~ 
ej ■ sgn [cos ft]. Turning to the noise, we may again de- 
compose Si(u>) into a frequency-independent part and a 
term arising from correlations between x(t) and m(t): 

Si(u) = £ diSi{i,u) = Sf + AS7(w). For the 
frequency- independent contribution Sf, we find: 



e V 1 , i\i 
— rl+{r'f 

7T 



n 



k B T eff (t)(x*(t)} - ^(Ax ) 2 
Zir 



(13)" 

where the bar indicates a time-average. The first term 
is the standard result for the shot noise of an AC-biased 
junction 21 . The second term indicates that the time- 
dependent motion of (x 2 (t)) (calculated from Eq. (0) 
can make a frequency-independent contribution to the 
noise. For v 3> Q, (x 2 ) responds only weakly to the 
time-dependence of T e ff(t), whereas for v ~ fi, the re- 
sponse becomes appreciable and 180 degrees out-of-phase 
with V(t). If in addition 70 <C 7, one finds a resulting 
suppression of the oscillator's contribution to Sf; this is 
shown in Fig. 2. Small resonances also occur when £1 is a 
multiple of v. Note that the oscillator modification of Sf 
is not captured by the classical picture of a fluctuating 
conductance. 

Finally, the frequency-dependent contribution ASi(u>) 
to the noise, which arises from correlations between x(t) 
and m(t), takes the simple form: 



AS7H = jYl As i(> y±l 
± 









H(S)] 


DC 





(14) 



where the omitted terms correspond to "quantum cor- 
rections" of the sort previously discussed. Without these 
terms, Eq. i|14|) is precisely the answer expected for a 
fluctuating classical conductance- one needs to simply 
shift the noise in the DC case up to the frequency v. In 
contrast, the quantum corrections to ASi for AC bias 
are not simply given by shifting the corresponding terms 
found for DC bias- one finds that the quantum correc- 
tions are larger in the AC case by a factor of 4/tt. The 
effect of ASi(uj) on the full noise is shown in Fig. 3. 

In conclusion, we have presented a fully quantum me- 
chanical calculation of the frequency-dependent current 
noise of a tunnel junction displacement detector, for both 
the cases of DC and AC voltage bias. The oscillator can 
lead to large effects in the shot noise, even if the coupling 
to the detector is weak; moreover, these effects cannot be 
completely described using a classical picture of a fluctu- 
ating junction conductance. We thank Konrad Lchncrt 
and Florian Marquardt for useful discussions; this work 
was supported by the NSF under grant nsf-itr 0325580 
and by the W. M. Keck Foundation. 
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